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Problem  Statement 


The  properties  of  several  discretizations  of  nonlinear  boun¬ 
dary  value  problems  in  ordinary  differential  equations  were  stu¬ 
died.  Special  attention  was  given  to  the  performance  of  these  me¬ 
thods  on  problems  arising  in  bifurcation  theory  and  in  optimal  con¬ 
trol.  The  methods  considered  were  primarily  collocation  with  piece- 
wise  polynomial  functions. 

Summary  of  the  Most  Important  Results 

Methods  introduced  by  Doedel  (SINUM  (1979))  were  generalized 
and  studied  on  general  meshes.  They  were  shown  [ l]  ,  [7]  ,  [llj  to 
be  unconditionally  stable  without  a  quasi- uniform  assumption.  Thus, 
these  methods,  which  can  be  implemented  very  efficiently  [  1  ]  could 
potentially  be  used  in  combination  with  a  mesh  selection  strategy. 

This  has  not  been  investigated  yet.  The  relation  between  the  num¬ 
ber  of  solutions  of  the  discrete  problem  and  the  number  of  solutions 
of  the  continuous  problem  was  also  investigated.  Examples  were  de¬ 
veloped  to  show  that  these  need  not  be  the  same  in  general,  not 
even  asymptotically  as  the  mesh  size  goes  to  zero.  Only  under  very 
restrictive  assumptions  on  the  differential  equation  can  these  mul¬ 
tiplicities  be  guaranteed  to  agree  for  sufficiently  small  mesh  sizes 
[7]  . 

Collocation  at  Gauss  points  was  shown  to  be  a  high  order  discre¬ 
tization  of  certain  unconstrained  optimal  control  problems.  Best 
possible  convergence  rates  and  superconvergence  results  were  estab¬ 
lished  in  [7]  .  The  use  of  these  methods  with  an  adaptive  mesh  se¬ 
lection  strategy  was  considered  in  [3]  .  These  results  are  promising, 
but  unfortunately  incomplete. 


Much  of  the  effort  under  this  grant  was  spent  analyzing  com¬ 
putational  methods  for  parameter  dependent  two-point  boundary  value 
problems.  The  effects  on  the  accuracy  of  collocation  and  certain 
finite  difference  discretizations  in  the  neighborhood  of  bifurcation 
points  were  studied  in  [s],  [s]  ,  [9]  and  [lo].  Included  here  were 
both  Hopf  and  simple  bifurcation.  The  behavior  of  Newton's  method 
as  a  computational  scheme  to  solve  the  algebraic  equations  generated 
by  the  discretizations  was  studied  in  C  4  3  and  is  now  well  understood. 
One  important  result  obtained  was  the  derivation  and  computational 
justification  C9]  of  a  reparameterization  of  periodic  problems  to 
allow  for  the  computation  of  branches  around  turning  points. 

Consider  the  system  du/dt  +  SF(a,u)  =  0  where  we  seek  211-periodic 
solutions.  To  find  a  neighboring  solution  (S,a,u)  a  distance  r  from 
any  given  solution,  add  the  condition 

61  (s“so)2  +  6 2  (a_a0)2  +  6 3  llu“vli  88  r2 

where  the  6^  are  fixed  weights  and  |  |*  | j  is  the  L  -norm.  The  solution 
is  unique  up  to  translation  and  can  be  anchored  by  requiring  (u(O)-v(O) 
to  be  orthogonal  to  v'(0). 
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